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Q^, Sharp fronts observed by Chandra satellite between dense cool cluster cores moving 

■^ I with near-sonic velocity through the hotter intergalactic gas, require strong suppres- 

Q^ I sion of thermal conductivity across the boundary. This may be due to magnetic fields 

tangential to the contact surface separating the two plasma components. We point out 
^ ■ that a super-Alfvenic motion of a plasma cloud (a core of a merging galaxy) through 

a weakly magnetized intercluster medium leads to " magnetic draping" , formation of a 
thin, strongly magnetized boundary layer with a tangential magnetic field. For super- 
^T , sonic cloud motion, Ms > 1, magnetic field inside the layer reaches near-equipartition 

\^ \ values with thermal pressure. Typical thickness of the layer is ~ L/M^ <^ L, where L 

is the size of the obstacle (plasma cloud) moving with Alfven Mach number M^ S> 1. 
To a various degree, magnetic draping occurs both for sub- and supersonic flows, ran- 
dom and ordered magnetic fields and it does not require plasma compressibility. The 
^ • strongly magnetized layer will thermally isolate the two media and may contribute to 

c/3 . the Kelvin-Helmholtz stability of the interface. Similar effects occur for radio bubbles, 

quasi-spherical expanding cavities blown up by AGN jets; in this case the thickness of 
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the external magnetized layer is smaller, ~ L/M\ <C L. 



1. Introduction 

Chandra observations of intergalactic medium (IGM) in clusters of galaxies often show sharp 
discontinuities in gas density, that separate dense cool gas moving with near-sonic velocity through 
the hotter gas (Markevitch et al. 2000; Vikhlinin et al. 2001). These fronts result from cluster 
mergers, when a cold subcluster core moves through a hot IGM at transonic velocities. The front 
is very sharp. For example, in the case of Abell A3667 the thickness of the front may be as thin as 
~ 5 kpc, (Vikhlinin et al. 2001). 

This sharpness is surprising since heat conduction together with development of the Kelvin- 
Helmholtz instability from tangential motion of gas are expected to result in much broader tran- 
sitions. Subtle geometrical and kinetic effects may stabilize the contact (Churazov & Inogamov 
2004) but this does not alleviate the problem of heat conduction. For typical cluster parameters, 



the Spitzer mean free path in unmagnetized plasma is ~ 30 kpc, much larger than the thickness of 
the transition layer. Thermal conductivity in unmagnetized gas should lead to cloud evaporation 
on time scale of ~ 10^ years (Ettori & Fabian 2000). Suppression of conductivity due to saturation 
of heat flux, when electron mean free path is comparable to the scale length of the temperature 
gradient (Cowie & McKee 1977), is not sufficient. This prompted Ettori &; Fabian (2000) to argue 
that classical conductivity has to be reduced by a factor of between 250 and 2500. 

Magnetic field is the prime suspect in reducing the conductivity across the front, since Larmor 
radius, typically in thousands of kilometers, is many orders of magnitude smaller than both electron 
mean free path and the scale of the temperature gradient. It was suggested that fields turbulent on 
small scales may do the job, reducing Spitzer conductivity by a large factor (Chandran & Cowley 
1998). On the other hand, Narayan & Medvedev (2001) argued that that if the scale of magnetic 
turbulence is smaller than the mean free path, and if the fluctuation spectrum extends over several 
decades in wave-vector, thermal conductivity is strongly enhanced almost up to the Spitzer value. 
A possible caveat in this argument in application to cold fronts is that it assumes that magnetic 
fields in the two media do intertwine. Vikhlinin et al. (2001) suggested that cross-boundary drift, 
heat conduction and KH instabilities are all suppressed by a large, > lOfiG, magnetic field along the 
contact surface. Such large, near-equipartition (with thermal pressure) magnetic fields must be local 
to the front in view of the estimates of sub-microgauss fields by Faraday rotation measurements in 
radio sources seen through the IGM (e.g. Kim et al. 1991), (see reviews by Carilli &: Taylor 2002; 
Enfilin et al. 2005; Govoni 2006). Vikhlinin et al. (2001) suggested that a magnetic field parallel 
to the contact surface arises due to shearing of an initially turbulent magnetic field . 

In this letter we discuss a straightforward dynamical effect that leads to formation of a narrow 
layer of tangential, near-equipartition magnetic field at the contact discontinuity. Such region forms 
regardless of however small the magnetic field is in the bulk of a cluster. This effect is well-known 
in space physics, and sometimes is called magnetic draping or magnetic barrier. It is mostly 
pronounced for the interaction of the Solar wind with Venus and Mars, planets that do not have 
their own magnetic field . It is also related to the so-called plasma depletion layer, where the plasma 
density is depressed with respect to values in the rest of the magnetosheath (Zwan & Wolf 1976; 
Paschmann et al. 1978). 



2. Properties of the transition layer 

2.1. Formation: divergence of a magnetic field on the contact in kinetic 

approximation 

To explain the effect of magnetic draping in its simplest form, consider an interaction of a 
cloud moving through a weakly magnetized medium parametrized by the ratio of total plasma to 
magnetic pressure (3 = Sirp/B'^ ^ 1. Since the magnetic field is weak, one may be prompted to 
neglect its dynamical effects on the flow completely. It turns out, that such approximation is not 



self-consistent and should break down close to the contact discontinuity separating two fluids. 

As a first step, which will be shown to be faulty, let us consider magnetic field in kinetic 
approximation, when dynamical effects of the magnetic field on the flow are neglected and the field 
lines are just advected with the flow satisfying froze-in condition. For simplicity, let us first assume 
that magnetic field in bulk is large scale and is directed along y axis, orthogonally to the direction of 
cloud motion in z direction (see Fig. 1). In addition, we assume that the cloud is axially symmetric 
around the stagnation line. The motion of magnetic field generates an inductive electric field, which 
under ideal approximation is E = — v x B. Combining induction equation 9jB = — curlE with 
continuity equation dtp + div pv = gives (Alfven 1942) 

(5. + (vV))^ = -|^=(^.v)v (1) 

P dt p \p ) 

On the other hand if (51 is an infinitesimal vector connecting two fluid elements along the flow line, 
its evolution with time is described by the same equation -^SV = (51 • V) v. This shows that quantity 
— evolves according to the length along flow line. Now, it is easy to see that as field line is wrapped 
around the contact, its length increases in proportion to the radius of the magnetosphere w at a 
given location, B/p ~ [B/ p)q{w/wq) where wq is the initial "impact parameter" of a given field 
line. On the contact Wq -^ and as a consequence the ratio {B/p) -^ cx) (Pudovkin & Semenov 
1985). 

To understand this in a different way, note that Eq. (1) can be solved using Cauchy's integral 
[e.g. Stern 1966): 

B ^ Bq • yWr 

P Po 

where V'^^^r is a derivative of the coordinate of a fiow element r(ro,t) with respect to initial 
coordinate tq, Bq is magnetic field at infinity. In curvilinear coordinates with scale factors hi Eq. 
(2) becomes 
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(summation on j only) . For an axially symmetric fiow the scale factor h, corresponding to the 
cyclic variable (p is proportional to roo, initial radial cylindrical coordinate, which is equal to the 
distance from the symmetry axis. For toroidal component of magnetic field , Eq. (3) gives 

Ba, f Ba,\ ( w 
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where we used the fact that the azimuthal angle remains constant, 0o = 0, so that q^- = 1. Eq. 
(4) shows that B^ tends to infinity on the contact surface corresponding to vjq -^ 0. 

Divergence of the toroidal component of magnetic field can also be shown explicitly by rewriting 
(j) component of the Eq. (1) (in spherical or cylindrical coordinate system) as 

{dt + (v V)) ^ = (5) 
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implying that B^jwp remains constant along the flow lines and is equal to Bq^^/wqPq {e.g. , for a 
constant magnetic field at infinity along y direction Bq^^ = Bqsvii(J)). 

To illustrate the divergence of the other component of the magnetic field parallel to the contact 
(besides -B^) we will assume a particular form of the contact and a type of the fluid motion. As 
an example, consider a spherical body of radius Rq moving with velocity vq subsonically through 
an incompressible fluid. The velocity flow in the frame associated with the body at a point defined 
by unit vector n is given by v = i?Q(3n(nvo — vo))/(2r^) — vq where r is a distance to the center 
of the body {e.g. Landau & Lifshits 1975). We can find an equation for flow lines tu/tuo = 
1/y^1 — {Rq/t)^, which gives dO/dOo = {vj /ti]q){zq/ z) ~ 0/Oq (last relation assumes small 9 and 
^o)- Here ^ is a polar angle of a spherical system of coordinates aligned with the direction of the 
motion, Fig. 1. Thus, dO/dOo diverges on the contact zuq -^ (or ^o ~^ 0) implying that the 
component Bg of the magnetic field diverges on the contact. Explicit calculations of the magnetic 
field structure for this problem confirm this conclusion (Bernikov & Semenov 1980). Incompressible 
motion past an axially symmetric body of arbitrary shape can be obtained from the solution for a 
spherical body through conformal mapping of the boundaries. 

It is easy to see that the divergence of d6/d0Q at ^o ^ occurs more generally than just in the 
case of an incompressible flow considered above. It follows from the fact that near stagnation line 
radial (in cylindrical coordinates) velocity increases linearly with distance from the line, v^ oc w, 
which implies that dO /OOq ~ /6q. 



The above derivations simplify along the stagnation line tu = 0. Then the induction equation 
gives B^Vz = const = BqVq. Since at the stagnation point Vz = 0, it follows that B.^ -^ oo. 

So far we have assumed that there are no discontinuities (shocks) in the flow. In case of a 
supersonic motion the toroidal component of the magnetic field will experience a jump at the shock 
and, in addition, the flow lines will generally experience a bend. This will introduce correction 
factors to the above relations but will not remove the divergence of a magnetic field on the contact. 
The subsonic incompressible flow considered above to prove the divergence of the tangential to 
the contact component of a magnetic field should be a reasonable approximation near the contact 
surface far downstream of a possible forward shock. 

The derivations above show that magnetic field is dynamically amplified on the contact. This is 
generally independent of whether the motion is supersonic or subsonic and is independent of plasma 
compressibility, though details of magnetic field amplification will be somewhat different in these 
cases. In addition, since governing equations can be written in terms of total derivatives along 
flow lines and thus are independent on the global structure of a magnetic field , tangled fields are 
also subject to the same amplification of the component parallel to the contact. The amplification 
occurs due to longitudinal stretching of magnetic field lines. Thus, a magnetized boundary layer is 
created in which magnetic field may reach near-equipartition values (for supersonic bulk motion) 
regardless of its value in the bulk (it should be non-zero, though). This effect is called magnetic 
draping or magnetic barrier. 



Note also that the 2-D case is very different from 3-D case. For example, in 3-D, the induction 
equation for stationary flow gives 

(v V)(rsin6l£;^) = (6) 

so that the quantity wE^ is constant along the flow line, zuE^ = zuqEq^^. ^ This implies that 
toroidal electric field vanishes on the contact, wq = 0. Contrary to this, in 2-D the electric field 
is constant everywhere, including on the contact, E = {BiVn — BnVi)ex, where coordinate I is 
along the generators of the contact surface from the axis of the symmetry and n is orthogonal 
to it. Divergence of magnetic field on the contact then follows from the requirement that normal 
components of both the velocity and magnetic field should vanish, «„ = i?„ = 0, implying Bi ^ oo 
(cf. Parker 1973, model for planar flows). 

For super- Alfvenic motion. Ma > 1, the only case when a magnetized layer does not form near 
the contact is when a cloud moves along the external magnetic field . In this case both B^ = 
and E^ = everywhere. In an axially symmetric, incompressible, stationary flow, velocity and 
magnetic field may be expressed in terms of flux functions ^ and P, v = V^/zu, B = VP/zu where 
P = Afpw and A^p is a toroidal component of vector potential. Condition E^ = then gives {e.g. 
in spherical coordinates) 

doPdr^ - de^drP = (7) 

implying that P = /(^) and B^ = Vrf, Bg = vqJ' . In particular, if magnetic field at infinity 
is constant in space then /' =constant, so that two flux functions are linearly related, P oc ^ . 
Since for a spherical obstacle w ^ = and vq is finite on the contact, magnetic field remains finite 
everywhere. 

The divergence of magnetic field on the contact is, of course, a consequence of the kinetic 
approximation. In full MHD case magnetic field is expected to be amplified to the point when 
its dynamical influence cannot be neglected anymore. For sonic and supersonic motion this cor- 
responds to near-equipartition (with total pressure, not just those of relativistic particles) fields. 
The layer itself must be modeled using MHD, not fluid equations. Inside the layer non-isotropic 
magnetic pressure will break the axial symmetry of the flow even for an axially symmetric obstacle. 
Stagnation point flow changes into stagnation line flow, where stream lines branch off not only at 
the stagnation point, but also at all points along the symmetry line (x = plane in Fig. 1). In 
case of a flow of plasma cloud, since position of the contact is determined by the pressure balance 
between two media, this will lead to deformation of the contact surface (by a fraction ~ 1/M^). 

A magnetized layer should also be associated with depletion of plasma. To see the reason 
for this, note that the total pressure, which is a sum of magnetic and gas pressures, should re- 
main approximately constant well behind a possible forward shock (neglecting possible temperature 
anisotropy, see below). Increasing magnetic pressure requires a decrease in gas pressures, which in a 



^Tangential component of electric field is continuous across tangential discontinuities, like shocks. 



polytropic fluid with adiabatic index 7 {e.g. for isentropic ideal gas) is accompanied by a decrease in 
density (approximately as p oc (/?/(! + (3))^''^ , decreasing with decreasing /3). Similarly, the plasma 
temperature should decrease as well, T oc (/?/(! + P))^'^~^'''^ . Thus, for a supersonic motion, the 
plasma density first increases at the forward shock and then decreases inside the magnetic barrier 
close to the contact (Wu 1992). 

Reality is a bit more complicated than this simple estimate since several competing plasma 
physics effects come into play. First, compression of plasma in a magnetic field creates anisotropy of 
plasma pressure due to conservation of adiabatic invariants (double adiabatic model of Chew et al. 
1956). When compressed normally to magnetic field lines, conservations of the adiabatic invariant 
p±/{Bp) = const and effective polytropic index j± = 2 would lead to density and transverse 
temperature T± increasing with a magnetic field. But strong plasma anisotropy leads to expulsion 
of plasma from high magnetic field regions due to mirror forces (magnetic bottling, reflection of 
particles from regions of high magnetic field , is the best known example of mirror forces). Second, 
onset of plasma instabilities may limit the pressure anisotropy {e.g. to T± ~ Tli). In case of planets 
interacting with the Solar wind, typically isotropic MHD model give a reasonable fit to magnetic 
field , pressure and temperature profiles (Denton & Lyon 1996; Pudovkin et al. 1999; Song & Russel 
2002). Behavior of the magnetized boundary layer in some aspects is opposite to the purely fluid 
case. For example, the former leads to density minimum on the contact while the later predicts 
density maximum (Lees 1964), (see Song & Russel 2002, for a recent review). 

Historically, magnetic draping effect was somewhat a surprise in modeling of Solar wind in- 
teraction with planets. It was expected that for small magnetization the flow may be computed 
from purely hydro dynamical equations, and a magnetic field may be added later using frozen-in 
condition (Spreiter et al. 1966). Using this prescription Alksne (1967) found that the magnetic 
field goes to infinity at the contact, especially strongly at plane containing the symmetry axis and 
magnetic field (y = plane in our notations). Zwan & Wolf (1976) (see also Southwood & Kivelson 
1995) calculated in details the properties of the magnetic barrier and predicted that it should be 
associated with depletion of plasma density, as magnetic field lines are stretched along the contact 
surface and plasma is allowed to flow along stretched magnetic field lines and, in addition, develop- 
ing of temperature anisotropy leads to magnetic mirror forces pushing plasma away from regions 
of high magnetic field . 

These theoretical ideas have been generally confirmed by direct satellite observations of the 
Terrestrial {e.g. Paschmann et al. 1978; Crooker 1979; KalUo et al. 1994; Wang et al. 2003), 
Cytherean (Biernatet al. 1999) and Martian (0ieroset et al. 2004) depletion layers. Overall, 
observations seem to be consistent with modern full MHD models {e.g. Kallio et al. 1998; Erkaev 
et al. 2000, and preceding references). Magnetic draping also occurs at the outer heliosphere and 
may be related to low frequency, ~ 3 kHz, radiation observed by the Voyager spacecraft (Cairns 
2004). In astrophysical setting these ideas have been touched upon theoretically by Kulsrud et al. 
(1965); Rosenau &: Frankenthal (1976) for the case of supernova expansion and by Lyutikov (2002) 
in the case of relativistic GRB outflows. A number of numerical experiments also saw formation 



of the magnetic barrier (Mac Low et al. 1994; Jones et al. 1996; Gregori et al. 2000; Asai et al. 
2004, 2005), producing results in agreement with the above theoretical estimates. In particular, in 
case of merging cluster cores, the work of Asai et al. (2005) clearly shows formation of magnetic 
barrier and its larger thickness in 2-D, as expected (see section 2.2). 



2.2. Thickness of the magnetic barrier 

For low magnetization, /? ^ 1, a flow may be separated in two regions: a bulk, where motion is 
nearly hydrodynamic, and a boundary layer, where effects of a magnetic field are important. Flows 
in the two regions must be matched at their boundaries. Thus, magnetic draping almost does not 
affect, for example, the location of the forward shock in front of the obstacle. 

Next we estimate a thickness of the transition layer. Both for subsonic and supersonic motion 
of the cloud through IGM, the motion near the critical point is strongly subsonic and can be 
considered incompressible. In this case the velocity field is v^ = — (3Vo/2)(tz7/L), Vz = ?>Vq{z/L) 
(Landau & Lifshits 1975). Then, along the stagnation line tz7 = magnetic field evolves according 
to 

B^ + 2zdzB^ = (8) 

giving B-^ oc Xj \fz. To estimate the thickness of the magnetized layer, note that inside the layer 
magnetic pressure becomes of the order of ram pressure, B^ ~ STrpf q (at this point magnetic forces 
will strongly affect the plasma flow). If the typical size of the plasma cloud is L and M^ = ^oo/^A,inf 
is Alfven Mach number defined in terms of Alfven velocity Wyi.inf a-t infinity, then, using Eq (8) 
we find 
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The value of plasma (3 inside the magnetized sublayer is /3j„ ~ (1 + Mg))/Mg, where Mg = Vao/cg 
is sonic Mach number at infinity. Thus, for supersonic motion, Mg > 1, a near-equipartition layer, 
/3 ~ 1, forms (see also Zwan & Wolf 1976). In case of subsonic motion plasma (3 inside the sublayer 
is much smaller than in the bulk, /3i„//3 ~ 1/M\ <C 1 (so that plasma is more strongly magnetized). 
Both for subsonic and supersonic motion, the thickness of the magnetized boundary layer is much 
smaller that the size of the cloud (for Ma > 1). Note, the dimensionality of the problem is an 
important issue. Repeating the above estimates for a 2-D flow, the thickness of the magnetized 
layer is Ar/L ~ 1/Myi (Erkaev et al. 1995), much larger than in the 3-D case. 

The time it takes to form the layer may be estimated from a condition that a swept-up magnetic 
flux is of the order of the magnetic flux through the layer. A layer forms rather quickly, after the 
cloud traversed a length / ~ L/Ma < L. An exception is when the cloud moves very slowly, 
sub-Alfvenically Ma < 1, in which case no magnetized layer forms anyway. 



3. Application to cluster cold fronts 

A well-studied case is Abell 3667, which we use below for numerical estimates (Vikhlinin & 
Markevitch 2002). In this case the density of the hot component is n/j ~ 8 x 10""^ cm ~^, its 
temperature T/j ~ 8 keV, front velocity ~ 1500 km s^^, corresponding to sonic Mach number 
Ms ~ 1. The Spitzer Coulomb free path is Ag ~ 30 kpc (Ettori & Fabian 2000). If an average 
magnetic field in the cluster is ~ 1/uG, the electron Larmor radius is only ~ 4 x 10^ cm. Thus, the 
Larmor radius is many orders of magnitude smaller than the collision length, so that even in case 
of highest possible cross- field diffusion (Bohm diffusion), cross- field conductivity is virtually zero. 

Assuming B = l/xG, the plasma /3 parameter in the bulk is then (5 ~ 400, Alfven Mach number 
is Ma ~ y'7/3/2M<j ~ 20, so that for the core of L ~ 500/cpc the thickness of the strongly magnetized 
boundary layer is ~ L/M\ ~ 1.25kpc. This is somewhat smaller than Chandra resolution for a 
typical cluster, but since the thickness depends strongly on the assumed IGM magnetic field , oc B^, 
it is possible that in some cases the layer may be resolved. 

If the IGM is permeated by large scale magnetic field , the magnetized layer provides a contri- 
bution to rotation measure of the order of 

RM . 15.ad/,.^ m" f^V (^^f (^Y" (^-) (10) 
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This is a relatively small value. In addition, since only the component of magnetic field along the 
line of sight contributes to the rotation measure, this estimate is subject to strong geometrical 
variations over the contact. 



4. Radio bubbles 

Magnetic draping should also be important for stabilization of rising radio bubbles. X-ray 
emission voids of up to 30 kpc in size, against Rayleigh- Taylor and Kelvin-Helmholtz instabilities 
and, similar to merging cores, suppression of thermal conductivity across the front. For example, 
(Robinson et al. 2004) (see also (Churazov et al. 2001; Briiggen &; Kaiser 2002; Jones &: De Young 
2005)) showed that in the absence of magnetic field, bubbles are disrupted by hydrodynamic 
instabilities and effects of thermal conduction. In order to stabilize the contact, locally magnetic 
field should be of the order of equipartition field, while the field in bulk is much smaller. The effect 
of magnetic draping provides exactly what is needed for stability and suppression of conductivity: 
near-equipartition magnetic field tangential to the contact (in case of nearly sonic or supersonic 
expansion). Hydrodynamically, the problem of expanding cavity blown by an AGN jet is similar 
to the one considered by Kulsrud et al. (1965) of a supernova remnant expanding into ISM (see 
also De Young 2003). 

Consider an impermeable sphere expanding radially with radius R{t) into incompressible 
medium permeated by constant magnetic field i?o- From the incompressibility condition we find 
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Vr = R' {t){R{t) / r)"^ . Inductive electric field is in </> direction and the induction equation gives 

dr{B0Vrr) 



dtBe 



r 



d,Br = _M!i!i^ (11) 

r sin^ 
Assuming that Eg depends on self-similar variable ^ = r/R{t) > 1, equation for Bg becomes 

B'eif - 1) + ^ = (12) 

which gives 

sing 
^^ - (1_ ^-3)1/3^0 

Br = -cos9{l-C^)^^^Bo (13) 

This implies that tangential component of magnetic field diverges on the contact ^ = 1 as oc 
(3(.^ — 1))~^'^. Using Eq. (13) we can find equation for field lines, (^^ — l)sin0 = zuo/R{t), 
see Fig. 2. Similarly to the case of translational motion, amplification of magnetic field on the 
contact will lead to formation of a magnetized boundary layer which thickness now is smaller by 
a factor Ma- Ar/R ~ M^' {cf. Kulsrud et al. 1965). MHD simulations generally confirm this 
picture: Robinson et al. (2004) and Jones &: De Young (2005) find that modest IGM magnetic 
fields can suppress thermal conductivity and stabilize the rising bubbles against disruption by fluid 
instabilities. 

Presence of a large scale magnetic field inside a bubble, a leftover from AGN pumping, may 
ease constraints on stability. An internal magnetic field may also be needed to compensate for 
high external pressure since bubbles appear to be in pressure equilibrium, but the absence of X-ray 
emission argues for a lower temperature in comparison with external IGM. An additional stabilizing 
effect may come from a pile-up of magnetic field inside the contact discontinuity. Typically, the 
entropy of gas in the cores is lower than at the outskirts, resulting in motion of gas between the 
core and the contact. This will similarly create a magnetized layer on the inside of the contact. 



5. Conclusion 

We point out that presence of a dynamically unimportant magnetic field in the bulk of the 
IGM leads to formation of a strongly magnetized boundary layer which undoubtedly affects the 
mechanical and thermodynamical coupling between plasmas of the cold merging cores or rising 
AGN blown bubbles on the one side and hot IGM plasma on the other side. To a different degree 
this occurs both around sub-sonically and super-sonically moving flows, for random and large-scale 
magnetic flelds and does not require plasma compressibility. Primarily, magnetic draping leads to 
strong suppression of thermal conductivity across the contact and may explain observed narrowness 
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of the transition layers. A boundary layer of near-equipartition magnetic field (for supersonic motion 
Ms ^ 1) is also expected to stabilize the KH instability of the contact, especially close to the critical 
point of the flow (in case of merging cores). 

Thus, even weak bulk magnetization strongly affects the interaction of a flow with an obstacle. 
This runs contrary to expectations that a weak magnetic field in the bulk should not affect much 
overall dynamics of a merging cluster. For example, Heinz et al. (2003) (see also Nagai & Kravtsov 
2003; Bialekei al. 2002) carried out hydrodynamical simulations of the interaction of cold subcluster 
plasma and hot ambient matter neglecting heat conduction and magnetic fields. When Asai et al. 
(2004) repeated simulations of Heinz et al. (2003), they found that when the Spitzer conductivity 
is adopted, the subcluster evaporates rapidly and the cold front is not formed. 

The main prediction of the model is that magnetic field may reach near-equipartition values 
and be directed along the contact separating two fluids. This is best tested with high resolution 
polarization radio observation. If high-energy non-thermal particles are accelerated locally and 
produce synchrotron emission, the direction of the magnetic field may be then inferred from linear 
polarization. For example, radio observations of NGC 4522, a prominent galaxy in the Virgo cluster, 
indeed show magnetic fields along the front (Vollmer et al. 2004). Parallel magnetic fields are also 
observed in the case of peripheral features seen in Abell 2256 (Rottgering et al. 1994; Clarke 
2001). In addition, a magnetized layer may contribute to the rotation measure of background radio 
sources. Naturally, in both cases we can estimate only average values of magnetic field along the 
line of sight and it's not easy to single out contributions from a narrow, strongly magnetized layer. 

A number of observations may already be interpreted in the framework of the model: enhanced 
RMs are often seen at the edges of structural features in clusters. E.g., Carilli et al. (1988) see 
enhanced RMs at the edges of hot spots in Cygnus A (Carilli et al. 1988), Taylor et al. (1992) 
find an enhanced RM at the edge of one of the hot spots in 3C 194. Though our estimates in 
application to Abell 3667 show that variations of RM across the front are typically small (and even 
smaller for expanding bubble due to smaller thickness in that case. Section 4), total RM is a strong 
function of assumed magnetic field in the bulk, ~ B^. The Expanded Very Large Array (EVLA) 
will be most instrumental in mapping magnetic fields in clusters of galaxies. 

Another prediction of the model is that due to the depletion of plasma from the magnetized 
sheath we may observe a narrow layer ~ Ikpc thick of suppressed X-ray emission on the outside 
edge of cold fronts. Naturally, observations needed to test this are very challenging. 

Dynamic amplification of an external magnetic field on both sides of the contact creates con- 
ditions favorable for reconnection between external and internal magnetic fields (e. 5. Pudovkin 
et al. 2002). As a result of reconnection, normal components of velocity at the contact surface 
may be non-zero, determined either by stripping or by physics of a reconnection layer. Reconnec- 
tion is known to be an efficient site of particle acceleration which may produce observable radio 
and non-thermal X-ray signals (c/. a model of non-thermal emission in young supernova remnants 
of Lyutikov & Polh 2004). Wrapping of magnetic field lines may create conditions favorable for 
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reconnection in the wake of the core, similar to reconnection in the Earth magneto-tail (Galeev 
1979). 

Numerical simulations required to correctly describe the structure of the magnetic barrier are 
bound to be complicated. First, they must be done in full 3-D using MHD approximation: fluid 
models lead to qualitatively different results, while 2-D MHD simulations grossly overestimate the 
thickness of the magnetic barrier (by a factor Ma ^ 1 if compared with the 3-D case). In addition, 
since it is expected that plasma inside the magnetic barrier becomes anisotropic, simulations should 
account for possible temperature anisotropy, e.g. within the framework of Chew-Goldberger-Low 
(Chew et al. 1956) theory. Finally, kinetic effects, like self-limiting temperature anisotropy, may 
be important as well. 

Finally we note that the effect of magnetic draping may be important in other astrophysical 
applications, in particular in case of suppressed conductivity between dense cold HI ISM clouds and 
a lower density warmer medium {e.g. review by Bregman 2004). Supernova shocks and outflows 
from OB associations induce a large scale motion in the warm medium which will lead to magnetic 
blanketing of cold clouds. This effect should be especially important for High- Velocity Clouds. 

I would like to thank Andrei Kravtsov, Jean Eilek, Maxim Markevitch, Christoph Pfrommer, 
Dmitry Uzdensky and Alexey Vikhlinin for comments and discussions. 
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Fig. 1. — Draping of magnetic field lines around an axially symmetric obstacle. Bq is magnetic 
field at infinity along y direction, a cloud is moving along z direction, dashed line is a flow line of a 
fluid element, wq is the initial distance from the symmetry axis, is a polar angle between a unit 
radius- vector from the origin n to the current position of a fluid element and z axis. 
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Fig. 2. — Magnetic field lines in the z — y plane for subsonic incompressible expansion of a sphere 
into a medium with a constant magnetic field along z axis at infinity, kinematic approximation. 
Note the strong compression of field lines on the contact r = R{t). 



